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Abstract 


This paper outlines some of the major unique features of rubber 
composite structures. The features covered are those related to 
the material properties, but the analytical features are also 
briefly discussed. It is essential to recognize these features 
at the planning stage of any long-range analytical, experimental, 
or application program. 

The development of a general and comprehensive program which fully 
accounts for all the important characteristics of tires, under 
all the relevant modes of operation, may present a prohibitively 
expensive and impractical task at the near future. There is therefore 
a need to develop application methodologies which can utilize 
the less general models, beyond their theoretical limitations and 
yet with reasonable reliability, by proper mix of analytical, 
experimental, and testing activities. 
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INTRODUCTION 


The textile-cord-reinforced-rubber composites used in various 
industrial products differ in many respects from the classical 
rigid composites. There are a great number of papers, books, and 
professional journals devoted to rigid composites. The literature 
on rubber composites, however, is very limited as compared to 
classical composites. To understand the mechanics of such composites, 
it is essential to develop an in-depth understanding of the way in 
which the internal variables of the constituents participate and 
interact in responding to external agents, i.e., mechanical, 
thermal or other environmental forces. It is, however, important 
to first study the properties of each constituent before dealing 
with the material properties of rubber composites. 


CORDS 


The textile cord reinforcements are structural members which may be 
viewed as one dimensional in a microscopic sense. The cords 
generally consist of several yarns twisted together. The yarns 
also consist of numerous filament components organized in a 
geometrical array with the view towards enhancing certain target 
properties. The cord properties therefore depend on the properties 
of the filaments and the geometrical organization as well as the 
interfacial characteristics. The filament itself has also been known 
to possess an internal structure at the microscopic level. 

In view of the above consideration, the cord itself is a complex 
structure which may be studied at the microscopic level (figs. 1 
and 2) . The way in which the cord properties are related to 
filament properties, geometries (ref. 1) and other variables is 
the subject of textile mechanics (fig. 3) . Properties such as 
strength, stiffness, and fatigue characteristics can be con- 
trolled by internal variables when their relationships are well 
understood. 

From a higher scale of continuum mechanics we are, however, 
concerned only with the phenomenological properties of the cords as 
experimentally obtained without any attempt for relating such 



measured values to the micro-structure in the sense of textile 
mechanics. The properties of interest at this continuum level are 
therefore the effective cord properties. 

Such an approach enables us to bypass the complexity of the textile 
mechanics in our formulation of composite properties. The 
limitation is, however, that such measured properties serve as 
"averages" only and therefore the continuum elements should be at 
least in the order of cord diameters. 



Figure 1. Typical tire cord. 
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(1) RAYON (5) NYLON -66 

(2) NOMEX (6) STEEL 

(3) POLYESTER (?) GLASS 

(4) NYLON- 6 


TESTED AT 70°F.,65% R.H. 
INSTRON - 100% STRAIN RATE 


Figure 3. Stress-strain characteristics of typical tire cords (2). 
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The cord is the major load-bearing member of rubber composite 
structures and as such should provide strength and many other 
characteristics of interest. Some of the expected performance 
characteristics of, the tire, and for that matter any other 
structure, can be directly related to cord properties. Let us 
consider some of the tire performance characteristics which are 
affected by cord properties. A partial list is given in Table 1. 


Table 1 

TIRE PERFORMANCE AFFECTED BY CORDS 
o Burst Strength 
o Bruise Resistance 
o Tire Endurance (Separations) 
o Power Loss 
o Tread Wear 
o High-Speed Endurance 
o Tire Size and Shape 
o Groove Cracking 

o Flat Spotting and Non-Uniformity 
o Tire Cornering Force 
o Tire Spring Rate 
Noise 


o 


Tables 2 to 4 of reference (3) provide a list of tire 
requirements and the related cord requirements. Such relations 
should be viewed with caution and qualifications. 


Vehicle Trends 
High Performance 

Downsizing 

Front Wheel Drive 


Fuel Economy 


Table 2 
(REF. 3) 

RADIAL PASSENGER 

Tire Requirements 

Low Aspect Ratio 

(Increased Cornering 
Forces ) 

Downsizing - Monoply 

Improved Tread Wear 
Retreading 

Rolling Resistance 


Cord Requirements 
Modulus 

Lateral Stiffness 

Dimensional Stability 

Fatigue 

Modulus 

Tenacity 

Dimensional Stability 
Lateral Stiffness 

Aged Adhesion 
Thermal Stability 
Fatigue 

Tenacity 

Modulus 

Hysteresis 
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Table 3 
(REF. 3) 
RADIAL TRUCK 


Vehicle Trends 

Increased Drive 
Position Loads 

Improved Fuel 
Economy/New 
Transportation 
Act 


Tire Requirements 
High Deflection 


Reduced Weight 


Retread 


Tire Requirements 
Increased Durability 

Lighter Weight 
Low Profile 
Retreading 

Rolling Resistance 


Cord Requirements 
Tenacity 

Strength/Area Ratio 

Tenacity 
Fatigue 
Aged Adhesion 
Fatigue 

Thermal Stability 

Tenacity 

Modulus 

Hysteresis 


Table 4 
(REF. 3) 

RADIAL AIRCRAFT 


Cord Requirements 

Tenacity 

Modulus 

Thermal Stability 
Tenacity 

Strength/Area Ratio 
Aged Adhesion 

Dimensional/Thermal Stability 
Fatigue 
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CORD VISCOELASTICITY 


The textile cords used in rubber-reinforced composites are often 
nonlinear viscoelastic. The viscoelastic deformations are associated 
with the loss of energy. The dissipated energy appears as heat and 
leads to temperature rise which in turn affects the material properties. 
The cord is therefore an important contributor to the energy loss in 
rubber-composite structures such as tires (4) . Figures 4 and 5 show 
some of the viscoelastic properties of polyester cords. 



Figure 4. Carpet plot of 1000/2 polyester real modulus. 
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SUMMARY ON CORDS DISCUSSION 


o Cords are structures and as such their effective properties are 
geometry and boundary condition dependent. 


o The effective cord properties are different when cords are 
embedded in the rubber matrix. 


o The calculation of the three-dimensional effective properties 
by analytical homogenization is very complex and impractical. 

o The properties of some cords such as nylon are highly temper- 
ature sensitive, particularly around the glass transition 
temperature. 


o The properties are different in tension and compression, 
o Cords are often nonlinear and viscoelastic. 
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RUBBER 


The physical properties of rubber compounds depend on various 
processing parameters and components. Here again, we are not 
concerned as how these properties are related to the molecular 
structure of the rubber and the physics of rubber vulcanization. 

The focus is on determination of the relevant properties by proper 
experiment and within the framework of continuum physics. The most 
important property is the elasticity of the rubber, which is 
distinctly different from other conventional materials. The most 
distinctive features of rubber elasticity are the deformability and 
the rapid recovery of the deformations when loads are removed. 
Rubber remains elastic at extension ratios of several hundred 
percent. Such elastic characteristics make the rubber unique in 
this respect. In fact, the major developments in the theory of 
large elastic deformations evolved around application to rubber 
elasticity. The rubber elasticity is an important subject in 
understanding finite-element analysis of such products. 


RUBBER ELASTICITY 


To understand rubber elasticity, we may first examine the 
thermodynamics of reversible processes. The first and the second 
laws of thermodynamics state that (5) 

dE = Tds + dw 

where E is the internal energy, T is the absolute temperature, S 
is the entropy and w is the work done on the system. Experimental 
work has shown that the rubber elasticity resides basically in the 
entropy term. The rubber elasticity therefore has an entirely 
different molecular origin than other elastic materials whose 
elasticities are primarily associated with the increase in internal 
energy through changes in molecular or atomic spaces. 
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Much work has been carried out to formulate rubber elasticity. 

One example of the molecular approach is one which considers the 
molecular chain length having Gaussian distribution. The elasticity 
parameters are calculated from such quantities as finite molecular 
length and molecular weight between crosslinks. The entropy change 
resulting from Gaussian theory leads to 

AS = -JNK + Xg + x 3 " 1) 

in which N is the number of network chains per unit volume 
and , ^2 < and *3 are the three principal extension ratios 

along the three mutually perpendicular axes of strain for pure 
homogeneous strains. 

The above equation provides a first order approximation to rubber 
elasticity but is not adequate over a large range of deformations. 
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PHENOMENOLOGICAL APPROACH 

An entirely different approach (6) is the phenomenological approach 
of continuum mechanics. In this approach, the existence of a strain 
energy function is postulated. It has been shown that such a strain 
energy function should depend on deformation gradients. The equa- 
tions for isotropic incompressible elastic media are as follows: 


W = W (I 1 , I 2 , Ig) 
I, = X? + X 2 + X? 


I = x 2 X 2 + x 2 x 2 + x 2 x 2 

a 2 12 A 1 A 3 a 2 a 3 


I 


3 


x 2 x 2 x 2 

a 1 a 2 a 3 


For the incompressible case 


W = W (I x , I 2 ) 


I 


3 


1 


In series form 


W = ee C.. (I - 3) 1 (I, - 3) i,j = 0, 1, 2, 

ij J ^ 
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An alternative form, suggested in (7) , is as follows 


W = W (X 1 ) + W (X 2 ) + W (X 3 ) 

The stress-strain relations in the principal directions (the 
stresses are force per deformed area) are 


a 


a 


1 


2 


=s 


a 


3 


3W 

31 , 


3W 

31 , 


3W 

31 , 


(2X 2 ) 2 + || 2X 2 (X 2 + X 2 ) + h(0) 

2 

( 2 X 2 ^ + "gY 2 X 2 ( X^ + Xg) + h(0) 

2 

(2X 2 ) + || 2X 2 (X 2 + X 2 ) + h(0) 

2 
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The one-dimensional stress— strain relation can be obtained from 
the three-dimensional relations as shown in figure 6. 


a = 2 (X 2 



,3W . 1__ 3W 

^ 3I 1 X 2 3I 2 



Figure 6. One-dimensional stress-strain relation. 
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It can be seen that the tangent compressive moduli increase 
significantly as the compressive strain increases. This feature 
realistically describes the rubber response in compression. In the 
finite-element treatment of rubber materials, the rubber is often 
modeled as linear, on the ground that the strains are small and 
that the linear constitutive law should then provide a reasonable 
approximation. This view is very disputable. It suffices to say 
that the compressive strains remain small due to stiffening of 
rubber in compression. This feature can only be handled properly 
by nonlinear constitutive laws. In the linear model of rubber, the 
material stiffness remains unchanged and therefore moderate 
compressive forces produce high compressive strains that are not 
seen in the actual structure. Such exceedingly large strains result 
in distortion of elements which quickly leads to severe numerical 
problems . 


A difficult problem, in dealing with rubber elasticity, arises from 
the incompressibility constraint condition. The incompressibility 
condition leads to certain simplifications in the exact analysis of 
the problem, basically because of reduction in the number of 
unknown parameters. Such is not, however, the case with the finite- 
element approach. In the variational formulation, for example, the 
Lagrange multiplier introduces an additional unknown scalar 
function into the finite-element formulation. This unknown should 
be accommodated at the element level. The procedure results in a 
significant increase in the number of total unknowns and also may 
result in an ill-conditioned stiffness matrix. A great deal of 
research has been carried out to find the most suitable element for 
handling incompressibility. The subject is still open but the fact 
remains that the incompressibility imposes an additional burden in 
numerical analysis. The incompressibility condition, aside from 
being inconvenient in the finite-element analysis, is an 
approximation for rubber-like materials. Such an approximation 
becomes increasingly less accurate as the percentage of carbon 
black increases in the rubber compound. The exact enforcement 
of incompressibility is therefore not actually needed. 
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NEARLY INCOMPRESSIBLE MATERIALS 


The strain energy function for nearly incompressible materials can 
be obtained by a series expansion of the strain energy function 
about (13-1) r and retaining the leading terms of up to second order 

in (i 3 _ 1 ), as follows 


w = w i (I 1 , I 2 ) + W 2 (I 1 , I 2 ) . (I 3 - 1) 

+ Wg (l ± , I 2 ) . (Ig - l) 2 + . . . . 

|lg - 1| « 1 

W = ll (I 1 - 3) 1 (I 2 - 3) j + W 2 - (Ig - 1) + W 3 (I 3 - l) 2 + . . . 


One special case of the above eguations is when the function W 2 
and W 3 are considered constants. The constants are described as 
follows : 


W 2 ( 1 ^ I 2 > = H, 


^3 ^ 1 ’ * 2 ^ ^2 


H i " - (C io + 2 C 01 ) 
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The near incompressibility can now be enforced by assigning large 
values to H. In fact, as H approaches infinity, ( 1 3 _ ) approaches 
zero so that the strain energy remains finite. The higher the 
value of H, the closer the incompressibility would be satisfied. 

H is referred to as the penalty number. In finite-element analysis, 
however, the large values of H can lead to overriding stiffness 
which results in numerical problems. The penalty method nevertheless 
permits the satisfaction of near incompressibility without increasing 
the number of unknowns in finite-element formulation. Figure 7 shows 
a set of typical properties for the rubber. 



Figure 7. Stress-strain relation, (8). 
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FRACTURE AND FATIGUE PROPERTIES 


The fatigue of the rubber has been the subject of many investigations 
in the past. The rubber fatigue is intimately related to the nature 
of rubber fracture and cut growth. The fracture mechanics approach 
for rubber was first adopted by Rivlin and Thomas (9) and Thomas (10) 
who promoted the concept of the tear energy in describing the cut 
growth mechanism. The tear energy approach has been applied to the 
study of the crack growth problem and to the description of fatigue 
behavior of rubber. We only consider the mechanically induced fatigue 
and this therefore excludes the fatigue caused by or resulting from 
non-mechanical sources such as aggressive environment. Some of the 
non-mechanical sources, such as ozone cracking in elastomers, may, 
however, be more damaging than mechanical sources. 

Busse's early results (11) on NR compounds, shown in Table 5, clearly 
demonstrate the unique feature of rubber fatigue. 

Table 5 

EFFECT OF STRAIN CYCLE 


STRAIN CYCLE 


LIFE, MINUTE 


5 - 60% 

Extension 

10 

10 - 60% 

•1 

25 

15 - 60% 

II 

90 

17 - 60% 

II 

150+ 

20 - 60% 

00 

00 


Some experimental fatigue data on various compounds are shown in 
Figures 8 to 11 taken from References (12-14). 
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Figure 8. Fatigue life of NR as a function of minimum strain 

range (from ref. 12). 


FATIGUE LIFE 



Figure 9. Imposed extension versus number of cycles to break 

( from ref . 12) . 
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I0 8 



•/.MINIMUM STRAIN 

Figure 10. Cycle of failure versus minimum strain 
for NR and SBR (from refs. 13 and 14). 



% MINIMUM STRAIN 

Figure 11. Cycle of failure versus minimum strain for 
various compositions (from refs. 13 and 14). 



SUMMARY ON RUBBER 


o Rubber remains elastic at extension ratios of several hundred. 


o Rubber is almost incompressible. 


o The rubber elasticity can be most conveniently expressed 
through the strain energy function. 


o Rubber has unique fatigue properties. Extensive experimental 
work needs to be carried out for fatigue characterization of 
various rubber compounds. 


o Appropriate failure theories are needed for interpretation of 
stress data from finite-element analysis. 

o Properties are temperature dependent. 
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CORD REINFORCED RUBBER COMPOSITE - SINGLE PLY 


In the macro-mechanics approach to composites, the actual 
heterogeneous medium, Figure 12, is replaced by an "equivalent" 
homogeneous medium. The equivalent medium has the same geometry as 
the heterogeneous medium but different material properties. The 
"effective" properties of the latter may belong to a different 
class of symmetry than those of the constituent materials. For 
instance, for isotropic constituents, the effective properties may 
be isotropic, orthotropic or anisotropic depending on the internal 
geometry of the composite. 

Let us consider a two-dimensional composite made of an isotropic 
matrix and reinforced by doubly periodic arrays of isotropic 
fibers as shown in Figure 13. A composite representation of the 
problem is a homogeneous medium with the same external geometry and 
boundary conditions as those of the heterogeneous case, but with 
orthotropic properties. The effective orthotropic properties are 
then defined for a "representative" element. It is very important 
to first understand what is meant by that. Hashin (15) defined and 
discussed this question and has drawn some analogies to the case of 
homogeneous continua. In such continua, the hypothesis is that the 
continua retain their properties even for infinitesimal elements. 

Due to the discrete nature and microstructure of homogeneous media, 
the infinitesimal element of continuum mechanics should be large 
compared to the scales of the microstructure. The infinitesimal 
element therefore exhibits properties which are some statistical 
averages of the microstructure properties. On the other hand, the 
infinitesimal elements should be very small compared to the 
dimensions of the continua. It then follows that the physical 
quantities at a point, such as moduli, stress and strain components, 
are in reality associated with averages over infinitesimal elements 
and not with a geometrical point. Due to the complex nature of the 
microstructure, the properties of the continuum infinitesimal 
elements are determined experimentally as opposed to calculating 
the averages from a microstructural theory. The same comments also 
apply to the "representative" element of composite theories, that 
is to say that a "representative" element is the infinitesimal 
element of composite materials. It therefore should be large 
compared to the dimensions of its material phases. A 
"representative" element, as defined, would retain and represent 
the properties of the composite continua and furthermore these 
properties would be insensitive to boundary conditions. 


233 





A distinction must be made between a "representative" element and 
unit elements or unit cells. The latter is defined as building 
blocks, so that the continua can be constructed by repeated use of 
such units. For example, consider the composite of Figure 13. 

There are a number of unit cell configurations which can equally 
serve as building blocks. Some of these possible choices are shown 
in Figure 14. Figure 15 shows the boundary effects on shear 
deformations. The average properties of these units, unlike 
those of a "representative" element, are highly boundary condition 
dependent (16). A "representative" element, however, consists of a 
large number of unit cells. 

From the definition of a "representative" element, it is apparent 
that the properties calculated or experimentally measured would be 
rather insensitive to boundary conditions. The problem, however, 
is that such a calculation would be a formidable task. Most 
published works dealing with a calculation of effective properties 
use a unit cell as the basis of their computations. Even for unit 
cells, the exact solutions can be obtained for only some simple 
geometries and simple boundary conditions. Another approach 
adopted _ is to determine the upper and the lower bounds for these 
properties through approximate solutions of energy formulations. 
These bounds are, however, far apart for composites with high 
cord-to-matrix stiffness ratios, such as rubber composites, and are 
therefore of little practical use. 

The immediate question is, therefore, how sensitive these properties 
are to boundary conditions specified over the surface of the 
unit cell, and, furthermore, how the element boundary conditions are 
influenced by global boundary conditions. These questions have not 
been fully investigated in the literature. It can, however, be 
stated that such sensitivities should depend on the relative 
stiffnesses of the cord and the matrix. As the stiffness of the 
matrix approaches that of the reinforcing materials, the 
effective properties should become less sensitive to boundary 
conditions and obviously independent of the boundary conditions in 
the limiting case of identical constituents' properties. These 
sensitivities are thus of greater concern in rubber composites than 
conventional rigid composites. The ratio of cord to rubber 
stiffness can exceed 30,000 in some composites which is far greater 
than those of rigid composites. 


235 










SUMMARY ON BASICS OF RUBBER COMPOSITES 


o The calculation of effective properties of rubber composites 
by homogenization is subject to more limitations than rigid 
composites. These limitations need to be established on a 
sound theoretical basis. 


o In experimental determination of the effective properties, 
the effects of size and the boundary conditions should be 
investigated . 


o The homogenized properties of a single ply cannot be easily 
calculated since all the 3-D homogenized cord properties are 
not generally available, with sufficient accuracy. Some of 
the effective composite properties, however, may prove to be 
somewhat insensitive to inaccuracies in cord properties. 


o A necessary characteristic of a composite material is 

statistical homogeneity. A representative element, used for 
calculation of the effective properties, should be large 
compared to typical phase regions. A representative element 
of a single ply does not have such a character istic in the 
thickness direction. The out of plane properties of a single 
ply, therefore, are subject to question. 


o The effective properties are different in tension and 
compression . 
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CHARACTERIZATION OF LAMINATED COMPOSITES 


We now consider a laminated structure / composed of N layers of 
cord-reinforced composite materials, as shown in Figure 16. 

Each layer of heterogeneous composite may then be modeled as 
homogeneous but orthotropic with respect to the proper local 
coordinate of each layer. The constitutive equations of the 
laminated composite, however, must be obtained with respect to 
a global coordinate system XYZ, as shown in Figure 16. The 
transformation relations may then be used for the appropriate 
layers to carry out the required transformation. The layers 
are numbered from top to bottom and no symmetry is assumed with 
respect to any axis. The stress and moment resultants for the 
laminated structure in terms of stresses are defined as follows: 



Figure 16. Geometry of typical laminated composites. 
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The displacement components and the resulting constitutive relations 
are 


u = u 0 {x, y) + zk x (x, y) 
v = v 0 (x, y) + :k y (x, y) 
w = u (jf, y) 


- 


~ 







too 

Ox 

'V, 


Au 

/I 12 0 

0 

^ 16 

*,, 

*12 

*16 



A 1 2 
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A 26 
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*26 

dr n 
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Q. 
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where the laminate stiffnesses are related to the layer stillness by: 


A ,= 


V q‘*dz Lj= 1,2,6 

Jo L^i 

k = 1 


Bij = 


k = 1 


Clj'zdz ij= 1 , 2,6 


*w = 


C { : k :'z 2 dz /,./= 1,2.6 


C ,k 'k l dz i.j = 4,5 No summation on i~j 


Jl) 
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SUMMARY ON LAMINATED COMPOSITES 


o Each ply may be homogenized, subject to the same limitations 
as those of a single ply. 

o Complete homogenization is not possible since the coupling 
between the forces and the moments can only be accounted for 
by preserving nonhomogeneity in the thickness direction. 

o Classical kinematics constraints, such as the Kirchhoff 
hypothesis, do not apply. 

o Plate/shell type stiffnesses are extremely hard to measure. 

o Constitutive relations are nonlinear due to angle change 
between the cords of adjacent plies. 
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ANALYTICAL AND COMPUTATIONAL FEATURES 
KINEMATICS 

The rubber composite structures undergo large strains as well as 
large rotations. The major kinematic features are listed in table 6. 


Table 6 

KINEMATIC FEATURES 


LARGE RUBBER STRAINS 

LARGE ROTATIONS 

CORD ANGLE CHANGES 

NEAR INCOMPRESSIBILITY CONDITION 


CONSTITUTIVE RELATIONS FOR COMPOSITES 


Single ply 

The single-ply composites can be considered as linear orthotropic 
when referred to axes of symmetry. For large displacements but 
moderately small strains, the rubber and the cords can be considered 
linear, but the strain-displacement gradient relations are nonlinear. 
The most appropriate form of constitutive equations, in such case, is 

fc S. . = fc C. fc e. , 
o in o ljkl o kl 
t 

where Q C ijkl are orthotropic properties referred to coordinates 
of initial material symmetry. The term o S ij are tlie com P onents of 
the second Piolla-Ki rchhof f stress tensor and are the compo- 

nents of the Lagrangian strain tensor. This form is invariant under 
a rigid body motion and therefore needs to be updated due to the cord 
angle change. This does not hold for constitutive equations 
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which utilize other stress and strain measures. The major drawback 
of modeling each ply separately is the increased size of the 
finite-element problem. 


Several plies 

It is often more convenient to lump several plies together in one 
element and therefore reduce the size of the problem. The preceding 
equation can be utilized but not longer remains unchanged due to the 
change in the angle between the minus and the plus ply. In such 
models the orthotropic properties of the combined plies continually 
change as functions of the cord angle. 


Composites with nonlinear constituents 

It is often necessary to account for the rubber nonlinearities in 
the finite-element analysis, even for small strains. The reason 
for such a need is that the typical compressive forces can produce 
large rubber strains if the rubber stiffening in not properly 
accounted for in the material modeling. These unrealistic large 
strains may lead to element distortion and eventual loss of 
numerical stability in the finite-element model. No rigorous 
nonlinear composite constitutive equations have appeared in the 
literature for rubber composites. This awaits further research 
in this field. 


Various modeling levels may be required depending on the nature 
of the problem. These models and their respective features are 
listed in Table 7. 
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Table 7 


MODELING FEATURES 


ELEMENT TYPE FOR COMPOSITES 


MEMBRANE ELEMENTS 


TWO- AND THREE-DIMENSIONAL 
ELEMENTS REPRESENTING: 


SINGLE PLY 


SEVERAL PLIES 




SHELL-TYPE ELEMENTS 


FEATURES 


REQUIRES 2-D MATERIAL PROPERTIES, 
INTERPLIES MODELED INDEPENDENTLY 
ONE ELEMENT PER PLY AND ONE ELEMENT 
PER INTERPLY IN THE THICKNESS 
DIRECTION, NO NEED FOR MATERIAL 
UPDATE DUE TO THE CORD ANGLE CHANGE 


REQUIRES 3-D MATERIAL PROPERTIES, 
ONE ELEMENT PER PLY IN THICKNESS 
DIRECTION, NO NEED FOR MATERIAL 
UPDATE DUE TO THE CORD ANGLE 
CHANGE 


REQUIRES 3-D MATERIAL PROPERTIES, 
LESS ELEMENTS, MATERIAL PROPERTIES 
MUST BE UPDATED DUE TO THE CORD 
ANGLE CHANGE, CANNOT MODEL THE 
COUPLING BETWEEN THE BENDING 
AND IN-PLANE FORCES 


INCREASED MATERIAL PROPERTY INPUT, 
LEAST NUMBER OF ELEMENTS, 

CAN MODEL THE COUPLING BETWEEN THE 
BENDING AND IN-PLANE FORCES, 
KINEMATICS CONSTRAINTS OF THE 
CORRESPONDING SHELL THEORY, MATERIAL 
PROPERTIES MUST BE UPDATED DUE TO 
THE CORD ANGLE CHANGE 
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SUMMARY OF MAJOR ANALYTICAL AND COMPUTATIONAL NEEDS AND FEATURES 
o Large elastic strains, 
o Large rotations. 

o Incompressibility condition, natural rubber can be considered 
as incompressible, but filled rubber exhibits compressibility. 

o Nonconservative loading. 

o Element type and aspect ratio. 

o Self-adopt i ve schemes for load increments and step size, 
o Contact algorithms for frictional loadings, 
o Finite-element formulation in rotating coordinate system, 
o Substructuring for localized analysis. 

o Stresses are very erratic at regions of sudden change in 
stiffness, such as cord-rubber interface. When calculated 
from the finite-element displacement method, proper smoothing 
algorithms must be developed for nonlinear problems. 

o Sensitivity analysis for uncertain input, material properties 
or other variables. 
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